The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

Phylogeography and demographic history of macaques, fascicularis species group, in East Asia: Inferred from multiple genomic markers.

Climate changes at larger scales have influenced dispersal and range shifts of many taxa in East Asia. The fascicularis species group of macaques is composed of four species and is widely distributed in Southeast and East Asia. However, its phylogeography and demographic histories are currently poorly understood. Herein, we assembled autosomal, mitogenome, and Y-chromosome data for 106 individuals, and combined them with 174 mtDNA dloop haplotypes of this species group, with particular focus on the demographic histories and dispersal routes of Macaca fuscata, M. cyclopis, and M. mulatta. The results showed: (1) three monophyletic clades for M. fuscata, M. cyclopis, and M. mulatta based on the multiple genomics analyses; (2) the disparate demographic trajectories of the three species after their split ∼1.0 Ma revealed that M. cyclopis and M. fuscata were derived from an ancestral M. mulatta population; (3) the speciation time of M. cyclopis was later than that of M. fuscata, and their divergence time occurred at the beginning of "Ryukyu Coral Sea Stage" (1.0-0.2 Ma) when the East China Sea land bridge was completely submerged by the sea level rose; and (4) the three parallel rivers (Nujiang, Lancangjiang, and Jinshajiang) of Southwestern China divided M. mulatta into Indian and Chinese genetic populations ∼200 kya. These results shed light on understanding not only the evolutionary history of the fascicularis species group but also the formation mechanism of faunal diversity in East Asia during the Pleistocene.

Hybridization and its impact on the ontogenetic allometry of skulls in macaques.

The role of hybridization in morphological diversification is a fundamental topic in evolutionary biology. However, despite the accumulated knowledge on adult hybrid variation, how hybridization affects ontogenetic allometry is less well understood. Here, we investigated the effects of hybridization on postnatal ontogenetic allometry in the skulls of a putative hybrid population of introduced Taiwanese macaques (Macaca cyclopis) and native Japanese macaques (Macaca fuscata). Genomic analyses indicated that the population consisted of individuals with varying degrees of admixture, formed by male migration from Japanese to Taiwanese macaques. For overall skull shape, ontogenetic trajectories were shifted by hybridization in a nearly additive manner, with moderate transgressive variation observed throughout development. In contrast, for the maxillary sinus (hollow space in the face), hybrids grew as fast as Taiwanese macaques, diverging from Japanese macaques, which showed slow growth. Consequently, adult hybrids showed a mosaic pattern, that is, the maxillary sinus is as large as that of Taiwanese macaques, while the overall skull shape is intermediate. Our findings suggest that the transgressive variation can be caused by prenatal shape modification and nonadditive inheritance on regional growth rates, highlighting the complex genetic and ontogenetic bases underlying hybridization-induced morphological diversification.

Prednisolone and rapamycin reduce the plasma cell gene signature and may improve AAV gene therapy in cynomolgus macaques.

Adeno-associated virus (AAV) vector gene therapy is a promising approach to treat rare genetic diseases; however, an ongoing challenge is how to best modulate host immunity to improve transduction efficiency and therapeutic outcomes. This report presents two studies characterizing multiple prophylactic immunosuppression regimens in male cynomolgus macaques receiving an AAVrh10 gene therapy vector expressing human coagulation factor VIII (hFVIII). In study 1, no immunosuppression was compared with prednisolone, rapamycin (or sirolimus), rapamycin and cyclosporin A in combination, and cyclosporin A and azathioprine in combination. Prednisolone alone demonstrated higher mean peripheral blood hFVIII expression; however, this was not sustained upon taper. Anti-capsid and anti-hFVIII antibody responses were robust, and vector genomes and transgene mRNA levels were similar to no immunosuppression at necropsy. Study 2 compared no immunosuppression with prednisolone alone or in combination with rapamycin or methotrexate. The prednisolone/rapamycin group demonstrated an increase in mean hFVIII expression and a mean delay in anti-capsid IgG development until after rapamycin taper. Additionally, a significant reduction in the plasma cell gene signature was observed with prednisolone/rapamycin, suggesting that rapamycin's tolerogenic effects may include plasma cell differentiation blockade. Immunosuppression with prednisolone and rapamycin in combination could improve therapeutic outcomes in AAV vector gene therapy.

Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Prevalent and Disseminated Recombinant and Wild-Type Adeno-Associated Virus Integration in Macaques and Humans.

Integration of naturally occurring adeno-associated viruses (AAV; wild-type AAV [wtAAV]) and those used in gene therapy (recombinant AAV [rAAV]) into host genomic DNA has been documented for over two decades. Results from mouse and dog studies have raised concerns of insertional mutagenesis and clonal expansion following AAV exposure, particularly in the context of gene therapy. This study aimed to characterize the genomic location, abundance, and expansion of wtAAV and rAAV integrations in macaque and human genomes. Using an unbiased, next-generation sequencing-based approach, we identified the genome-wide integration loci in tissue samples (primarily liver) in 168 nonhuman primates (NHPs) and 85 humans naïve to rAAV exposure and 86 NHPs treated with rAAV in preclinical studies. Our results suggest that rAAV and wtAAV integrations exhibit similar, broad distribution patterns across species, with a higher frequency in genomic regions highly vulnerable to DNA damage or close to highly transcribed genes. rAAV exhibited a higher abundance of unique integration loci, whereas wtAAV integration loci were associated with greater clonal expansion. This expansive and detailed characterization of AAV integration in NHPs and humans provides key translational insights, with important implications for the safety of rAAV as a gene therapy vector.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation.

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

Diet-related factors strongly shaped the gut microbiota of Japanese macaques.

Although knowledge of the functions of the gut microbiome has increased greatly over the past few decades, our understanding of the mechanisms governing its ecology and evolution remains obscure. While host genetic distance is a strong predictor of the gut microbiome in large-scale studies and captive settings, its influence has not always been evident at finer taxonomic scales, especially when considering among the recently diverged animals in natural settings. Comparing the gut microbiome of 19 populations of Japanese macaques Macaca fuscata across the Japanese archipelago, we assessed the relative roles of host genetic distance, geographic distance and dietary factors in influencing the macaque gut microbiome. Our results suggested that the macaques may maintain a core gut microbiome, while each population may have acquired some microbes from its specific habitat/diet. Diet-related factors such as season, forest, and reliance on anthropogenic foods played a stronger role in shaping the macaque gut microbiome. Among closely related mammalian hosts, host genetics may have limited effects on the gut microbiome since the hosts generally have smaller physiological differences. This study contributes to our understanding of the relative roles of host phylogeography and dietary factors in shaping the gut microbiome of closely related mammalian hosts.

Isolation of novel simian adenoviruses from macaques for development of a vector for human gene therapy and vaccines.

IMPORTANCE: Adenoviruses are widely used in gene therapy and vaccine delivery. Due to the high prevalence of human adenoviruses (HAdVs), the pre-existing immunity against HAdVs in humans is common, which limits the wide and repetitive use of HAdV vectors. In contrast, the pre-existing immunity against simian adenoviruses (SAdVs) is low in humans. Therefore, we performed epidemiological investigations of SAdVs in simians and found that the SAdV prevalence was as high as 33.9%. The whole-genome sequencing and sequence analysis showed SAdV diversity and possible cross species transmission. One isolate with low level of pre-existing neutralizing antibodies in humans was used to construct replication-deficient SAdV vectors with E4orf6 substitution and E1/E3 deletion. Interestingly, we found that the E3 region plays a critical role in its replication in human cells, but the absence of this region could be compensated for by the E4orf6 from HAdV-5 and the E1 expression intrinsic to HEK293 cells.

Comparative genomics reveals the hybrid origin of a macaque group.

Although species can arise through hybridization, compelling evidence for hybrid speciation has been reported only rarely in animals. Here, we present phylogenomic analyses on genomes from 12 macaque species and show that the fascicularis group originated from an ancient hybridization between the sinica and silenus groups ~3.45 to 3.56 million years ago. The X chromosomes and low-recombination regions exhibited equal contributions from each parental lineage, suggesting that they were less affected by subsequent backcrossing and hence could have played an important role in maintaining hybrid integrity. We identified many reproduction-associated genes that could have contributed to the development of the mixed sexual phenotypes characteristic of the fascicularis group. The phylogeny within the silenus group was also resolved, and functional experimentation confirmed that all extant Western silenus species are susceptible to HIV-1 infection. Our study provides novel insights into macaque evolution and reveals a hybrid speciation event that has occurred only very rarely in primates.

Genetic characterization of Macaca arctoides: A highlight of key genes and pathways.

When compared to the approximately 22 other macaque species, Macaca arctoides has many unique phenotypes. These traits fall into various phenotypic categories, including genitalia, coloration, mating, and olfactory traits. Here we used a previously identified whole genome set of 690 outlier genes to look for possible genetic explanations of these unique traits. Of these, 279 genes were annotated miRNAs, which are non-coding. Patterns within the remaining outliers in coding genes were investigated using GO (n = 370) and String (n = 383) analysis, which showed many interconnected immune-related genes. Further, we compared the outliers to candidate pathways associated with M. arcotides' unique phenotypes, revealing 10/690 outlier genes that overlapped these four pathways: hedgehog signaling, WNT signaling, olfactory, and melanogenesis. Of these, genes in all pathways except olfactory had higher FST values than the rest of the genes in the genome based on permutation tests. Overall, our results point to many genes each having a small impact on phenotype, working in tandem to cause large systemic changes. Additionally, these results may indicate pleiotropy. This seems to be especially true with the development and coloration of M. arctoides. Our results highlight that development, melanogenesis, immune function, and miRNAs may be heavily involved in M. arctoides' evolutionary history.

Systemic administration of novel engineered AAV capsids facilitates enhanced transgene expression in the macaque CNS.

BACKGROUND: Adeno-associated virus (AAV) vectors are a promising vehicle for noninvasive gene delivery to the central nervous system via intravenous infusion. However, naturally occurring serotypes have a limited ability to transduce the brain, and translating engineered capsids from mice to nonhuman primates has proved challenging. METHODS: In this study, we use an mRNA-based directed-evolution strategy in multiple strains of mice as well as a de novo selection in cynomolgus macaques to identify families of engineered vectors with increased potency in the brain and decreased tropism for the liver. FINDINGS: We compare the transgene expression capabilities of several engineered vectors and show that while some of our novel macaque-derived variants significantly outperform AAV9 in transducing the macaque brain following systemic administration, mouse-derived variants-both those identified in this study and those reported by other groups-universally do not. CONCLUSIONS: Together, the results of this work introduce a class of primate-derived engineered AAV capsids with increased therapeutic potential and highlight the critical need for using appropriate animal models to both identify and evaluate novel AAVs intended for delivery to the human central nervous system. FUNDING: This work was funded primarily through an anonymous philanthropic gift to the P.C.S. lab at the Broad Institute of MIT and Harvard and by a grant from the Howard Hughes Medical Institute to P.C.S.

A partial duplication of an X-linked gene exclusive of a primate lineage (Macaca).

Gene duplication plays a significant role in evolution. Paralogous gene copies may be lost due to the successive accumulation of deleterious mutations or remain active in the genome. In this work, a partial duplication of an X-linked region in the Macaca genus is identified and explored. Genomic comparisons reveal that the duplication encompasses the genes encoding ornithine transcarbamylase (OTC) and retinitis pigmentosa GTPase regulator (RPGR), spanning over 0.1 Mb on the chromosome 9 of Macaca. According to our analyses, the duplicated region of chromosome 9 involves partial coding sequences of both OTC and RPGR genes. Analyses of the selective pressures did not reveal significant differences in the ratio between nonsynonymous and synonymous mutations (w<1), suggesting that no selective pressures were acting in the evolutionary process. Reports for a biological role regarding some partial duplications exist in the literature, therefore, although being rare events, partial duplications of functionally important genes are worthy of study so that their impact can be explored.

Sequence Analysis of Novel Staphylococcus aureus Lineages from Wild and Captive Macaques.

Staphylococcus aureus is a widespread and common opportunistic bacterium that can colonise or infect humans as well as a wide range of animals. There are a few studies of both methicillin-susceptible S. aureus (MSSA) and methicillin-resistant S. aureus (MRSA) isolated from monkeys, apes, and lemurs, indicating a presence of a number of poorly or unknown lineages of the pathogen. In order to obtain insight into staphylococcal diversity, we sequenced strains from wild and captive individuals of three macaque species (Macaca mulatta, M. assamensis, and M. sylvanus) using Nanopore and Illumina technologies. These strains were previously identified by microarray as poorly or unknown strains. Isolates of novel lineages ST4168, ST7687, ST7688, ST7689, ST7690, ST7691, ST7692, ST7693, ST7694, ST7695, ST7745, ST7746, ST7747, ST7748, ST7749, ST7750, ST7751, ST7752, ST7753, and ST7754 were sequenced and characterised for the first time. In addition, isolates belonging to ST2990, a lineage also observed in humans, and ST3268, a MRSA strain already known from macaques, were also included into the study. Mobile genetic elements, genomic islands, and carriage of prophages were analysed. There was no evidence for novel host-specific virulence factors. However, a conspicuously high rate of carriage of a pathogenicity island harbouring edinB and etD2/etE as well as a higher number of repeat units within the gene sasG (encoding an adhesion factor) than in human isolates were observed. None of the strains harboured the genes encoding Panton-Valentine leukocidin. In conclusion, wildlife including macaques may harbour an unappreciated diversity of S. aureus lineages that may be of clinical relevance for humans, livestock, or for wildlife conservation, given the declining state of many wildlife populations.

Variants of the adeno-associated virus serotype 9 with enhanced penetration of the blood-brain barrier in rodents and primates.

The development of gene therapies for the treatment of diseases of the central nervous system has been hindered by the limited availability of adeno-associated viruses (AAVs) that efficiently traverse the blood-brain barrier (BBB). Here, we report the rational design of AAV9 variants displaying cell-penetrating peptides on the viral capsid and the identification of two variants, AAV.CPP.16 and AAV.CPP.21, with improved transduction efficiencies of cells of the central nervous system on systemic delivery (6- to 249-fold across 4 mouse strains and 5-fold in cynomolgus macaques, with respect to the AAV9 parent vector). We also show that the neurotropism of AAV.CPP.16 is retained in young and adult macaques, that this variant displays enhanced transcytosis at the BBB as well as increased efficiency of cellular transduction relative to AAV9, and that it can be used to deliver antitumour payloads in a mouse model of glioblastoma. AAV capsids that can efficiently penetrate the BBB will facilitate the clinical translation of gene therapies aimed at the central nervous system.

Cross-species evolution of a highly potent AAV variant for therapeutic gene transfer and genome editing.

Recombinant adeno-associated viral (AAV) vectors are a promising gene delivery platform, but ongoing clinical trials continue to highlight a relatively narrow therapeutic window. Effective clinical translation is confounded, at least in part, by differences in AAV biology across animal species. Here, we tackle this challenge by sequentially evolving AAV capsid libraries in mice, pigs and macaques. We discover a highly potent, cross-species compatible variant (AAV.cc47) that shows improved attributes benchmarked against AAV serotype 9 as evidenced by robust reporter and therapeutic gene expression, Cre recombination and CRISPR genome editing in normal and diseased mouse models. Enhanced transduction efficiency of AAV.cc47 vectors is further corroborated in macaques and pigs, providing a strong rationale for potential clinical translation into human gene therapies. We envision that ccAAV vectors may not only improve predictive modeling in preclinical studies, but also clinical translatability by broadening the therapeutic window of AAV based gene therapies.

Cognitive genomics of learning delay and low level of social performance monitoring in macaque.

Cognitive skills and the underlying neural architecture are under the influence of genetics. Cognitive genomics research explores the triadic relationship between genes, brain, and cognition, with its major strategy being genotype-driven. Here we show that an inverse strategy is feasible to identify novel candidate genes for particular neuro-cognitive phenotypes in macaques. Two monkeys, originally involved in separate psychological studies, exhibited learning delay and low levels of social performance monitoring. In one monkey, mirror neurons were fewer compared to controls and mu suppression was absent in the frontal cortex. The other monkey showed heightened visual responsiveness in both frontal cortex and dopamine-rich midbrain, with a lack of inter-areal synchronization. Exome analyses revealed that the two monkeys were most likely cousins and shared variants in MAP2, APOC1, and potentially HTR2C. This phenotype-driven strategy in cognitive genomics provides a useful means to clarify the genetic basis of phenotypic variation and develop macaque models of neuropsychiatric disorders.

Aging gut microbiota of wild macaques are equally diverse, less stable, but progressively personalized.

BACKGROUND: Pronounced heterogeneity of age trajectories has been identified as a hallmark of the gut microbiota in humans and has been explained by marked changes in lifestyle and health condition. Comparatively, age-related personalization of microbiota is understudied in natural systems limiting our comprehension of patterns observed in humans from ecological and evolutionary perspectives. RESULTS: Here, we tested age-related changes in the diversity, stability, and composition of the gut bacterial community using 16S rRNA gene sequencing with dense repeated sampling over three seasons in a cross-sectional age sample of adult female Assamese macaques (Macaca assamensis) living in their natural forest habitat. Gut bacterial composition exhibited a personal signature which became less stable as individuals aged. This lack of stability was not explained by differences in microbiota diversity but rather linked to an increase in the relative abundance of rare bacterial taxa. The lack of age-related changes in core taxa or convergence with age to a common state of the community hampered predicting gut bacterial composition of aged individuals. On the contrary, we found increasing personalization of the gut bacterial composition with age, indicating that composition in older individuals was increasingly divergent from the rest of the population. Reduced direct transmission of bacteria resulting from decreasing social activity may contribute to, but not be sufficient to explain, increasing personalization with age. CONCLUSIONS: Together, our results challenge the assumption of a constant microbiota through adult life in a wild primate. Within the limits of this study, the fact that increasing personalization of the aging microbiota is not restricted to humans suggests the underlying process to be evolved instead of provoked only by modern lifestyle of and health care for the elderly. Video abstract.

Single-cell transcriptomics of adult macaque hippocampus reveals neural precursor cell populations.

The extent to which neurogenesis occurs in adult primates remains controversial. In this study, using an optimized single-cell RNA sequencing pipeline, we profiled 207,785 cells from the adult macaque hippocampus and identified 34 cell populations comprising all major hippocampal cell types. Analysis of their gene expression, specification trajectories and gene regulatory networks revealed the presence of all key neurogenic precursor cell populations, including a heterogeneous pool of radial glia-like cells (RGLs), intermediate progenitor cells (IPCs) and neuroblasts. We identified HMGB2 as a novel IPC marker. Comparison with mouse single-cell transcriptomic data revealed differences in neurogenic processes between species. We confirmed that neurogenesis is recapitulated in ex vivo neurosphere cultures from adult primates, further supporting the existence of neural precursor cells (NPCs) that are able to proliferate and differentiate. Our large-scale dataset provides a comprehensive adult neurogenesis atlas for primates.